The influence of pairing on the nuclear matrix elements of the neutrinoless [3(3 decays 
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We study in this letter the neutrinoless double beta decay nuclear matrix elements (NME's) 
in the framework of the Interacting Shell Model. We analyze them in terms of the total angular 
momentum of the decaying neutron pair and as a function of the seniority truncations in the nuclear 
wave functions. This point of view turns out to be very adequate to gauge the accuracy of the 
NME's predicted by different nuclear models. In addition, it gives back the protagonist role in this 
process to the pairing interaction, the one which is responsible for the very existence of double beta 
decay emitters. We show that low seniority approximations, comparable to those implicit in the 
quasiparticle RPA in an spherical basis, tend to overestimate the NME's in several decays. 
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The discovery of the massive character of the neutri- 
nos in the recent measurement at Super-Kamiokande [l[ , 
SNO and KamLAND Q, has openeded a new era in 
the neutrino physics. However, these experiments are 
sensitive only to the mass differences between the three 
neutrino species. Their absolute mass scale and hierarchy 
are still unknown. In addition, we don't know either if 
the neutrinos are Dirac or Majorana particles. The dou- 
ble beta decay is the rarest nuclear weak process. It takes 
place between two even-even isobars, when the decay to 
the intermediate nucleus is energetically forbidden due to 
the pairing interaction, that shifts the even-even and the 
odd-odd mass parabolas in a given isobaric chain. The 
two-neutrino decay is just a second order process in the 
weak interaction. It conserves the lepton number and 
has been already observed in several nuclei. A second 
mode, the neutrinoless decay Ov [3(3 can only take place 
if the neutrino is a Majorana particle and demands an 
extension of the standard model of electroweak interac- 
tions, because it violates the lepton number conservation. 
Therefore, the observation of the double beta decay with- 
out emission of neutrinos will sign the Majorana charac- 
ter of the neutrino and will establish the absolute mass 
scale of the neutrinos, hence deciding their mass hierar- 
chy. 

The expression for the neutrinoless beta decay half-life, 
in the + — > + case, can be brought to the following form 
BS: 

[T$>(0+- > 0+)]- 1 = Go, (m<°"> f^l)) 2 (1) 



M<°"> = M$> - (g-) 2 - M<? v) (2) 

where (m„) is the effective neutrino mass (a linear combi- 
nation of the neutrino mass eigenvalues whose coefficients 



are the corresponding elements of the neutrino mixing 
matrix), and Go„ is the kinematic phase space factor [6j. 

The important point at this stage is that, once the 
neutrinoless double beta decay be detected, to transform 
the measured half life in an accurate value of the effec- 
tive neutrino mass would require a precise computation 
of the nuclear matrix elements (NME's) of the decay op- 
erators. This, in turn, demands a detailed description of 
the structure of the nuclei involved in the process. A crit- 
ical analysis of the available predictions for the NME's of 
the potential 0i>(3{3 emitters (only about one dozen) was 
made recently by Bahcall et al. pj. Their conclusion 
was rather pessimistic, owing to the large dispersion of 
the calculated values. In a subsequent paper, Rodin et 
al. [8j have shown that many of the quasi-particle RPA 
(QRPA) calculations taken into account in Bahcall's sur- 
vey were obsolete, and that, when these are not consid- 
ered, the spread of the calculated values is much smaller. 
The aim of this work is to go one step further and to 
propose a much narrower band of values for the NME's, 
based in the predictions of large scale applications of the 
Interacting Shell Model (ISM) and in the analysis of the 
QRPA results in terms of the pairing content of their 
solutions. 

The matrix elements Mgy f T can be calculated in the 
closure approximation, that is good to better than 90% 
due to the large average energy of the virtual neutrino 
(^100 MeV) [9(. For the Gamow- Teller channel it reads, 

= (0+\h(\fi - ma, ■ a 2 )(q Q-)|0+) (3) 

and similar expressions hold for the other matrix ele- 
ments. h(\fl — r^|) is the neutrino potential (~l/r) ob- 
tained from the neutrino propagator. Higher order con- 
tributions (hoc) to the nuclear current produce the tensor 
term and add extra contributions to the Gamow- Teller 
expression of Eq. ([3]) [1(3] • 

Generically, the two body decay operators can be writ- 
ten in the Fock space representation as: 
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where the indices i, j, k, I, run over the single particle 
orbits of the spherical nuclear mean field. Applying the 
techniques of ref. we can factorize the operators as 
follows 



(5) 



The operators Pj^ annihilate pairs of neutrons coupled 
to J v in the father nucleus and the operators Pj„ sub- 
stitute them by pairs of protons coupled to the same J rir . 
The overlap of the resulting state with the ground state 
of the grand daughter nucleus gives the J^-contribution 
to the NME. The -a priori complicated- internal struc- 
ture of these exchanged pairs is dictated by the double 
beta decay operators. 

The ISM calculations reported in this letter are car- 
ried out in the spirit of the previous shell model works 

H EE El El- For the A = 76 and A=82 cases we make 



full calculations in the valence space (lp3/ 2 , 0f 5 / 2 , lpi/2> 
Ogg/ 2 ) using a newly built effective interaction that, start- 
ing with a G-matrix [HI ] has its matrix elements fitted to 
a large set of experimental data. For the A=124, A=128, 
A=130, and A=136 emitters, we make full calculations 
in the valence space (0g 7/2 , ld 5/2 , ld 3/2 , 2s 1/2 , 0h 11/2 ) 
with another interaction obtained in a similar manner. 
These model spaces and interactions will be discussed in 
detail elsewhere [13] • The dimensions of the shell model 
bases reach in some cases O(10 10 ). The present calcu- 
lations adopt the closure approximation, and model the 
short range and finite size corrections as in We use 
ro=1.2 fm to make the matrix elements dimensionless and 
g J 4=1.25. The choice of g^=1.25 instead of the quenched 
value g J 4=l-0 needed for the pure Gamow- Teller pro- 
cesses in nuclei is consistent with the use of the closure 
approximation, in which the multipole decomposition of 
the decay plays no role at all. In addition, even without 
closure, the need of a quenched gA in the J w =l + channel 
of the Ov decay, that has no reason to resemble the pure 
Gamow- Teller operator of the 2is decay, is not guaran- 
teed. Hig her order contributions to the nuclear current 
(hoc) [l0( are explicitly included for the first time in the 
ISM context, leading to reductions of the NME's in the 
15% range. Our final predictions for M 1 - ^ are gathered 
in table HI 

Except for doubly magic 48 Ca, whose NME is 
severely quenched, our predictions cluster around a value 
M^ « 2. The upper bounds on the neutrino mass for 
a half life of 10 25 y, that incorporate the phase space fac- 
tors, show a mild preference for the potential emitters 



TABLE I: ISM predictions for the Qu double beta decay ma- 
trix elements, with and without higher order contributions to 
the nuclear current (hoc). The effective neutrino mass corre- 
sponds to T 1 = 10 25 y. 







M (0,/) (no hoc) 


M (0W 


(m„) 


48 Ca^ 


48 T j 


0.76 


0.59 


1.09 


76 Ge^ 


76 Se 


2.58 


2.22 


1.05 


82 Se^ 


82 Kr 


2.49 


2.11 


0.50 


124 Sn - 


i24 Te 


2.38 


2.02 


0.53 


12 8 Te _> 


128 Xe 


2.67 


2.26 


2.27 


130 Te -» 


130 Xe 


2.41 


2.04 


0.41 


136 Xe -> 


136 Ba 


2.00 


1.76 


0.47 



with A=82, 124, 130 and 136. The matrix elements are 
dominated by the Gamow- Teller contribution. The in- 
fluences of the restrictions in the valence space and of 
the choice of the effective interaction in the ISM NME's 
have been studied in In all, these should result in 

a 20% uncertainty of our predictions. Treating the short 
range correlations with a prescription softer than Jastrow 
might produce an increase of NME's, that we have not 
evaluated yet in the ISM context, but we do not expect 
it to go beyond 10-15%. 

In order to explore the structure of the 0vf3(3 two body 
transition operators, we have plotted in Figure [U the 
contributions to the 0z/ GT matrix element as a function 
of the J" of the decaying pair. The results are very sug- 
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FIG. 1: (color online) Contributions to the Gamow Teller 
matrix element of the 82 Se — * 82 Kr decay as a function of the 
of the transformed pair (no hoc). 



gestive, because the dominant contribution corresponds 
to the decay of J=0 pairs, whereas the contributions of 
the pairs with J>0 are either negligible or have opposite 
sign to the leading one. This behavior is common to all 
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the cases that we have studied, as can be seen in table 
III1 Notice that the cancellations are substantial. These 
features are also present in the QRPA calculations, in 
whose context they had been discussed in refs. [H, . 



TABLE II: J=0 vs J>0 pair contributions to the Gamow 
Teller matrix element (no hoc). 





M&? 


M^(J=0) 


M^(J>0) 




2.35 


5.59 


-3.24 


82 ge ^ 82 fo 


2.25 


5.32 


-3.07 


130 Te ^ 130 Xe 


2.12 


6.58 


-4.46 


136 Xe ^ 136 Ba 


1.77 


5.72 


-3.95 



To grasp better this mechanism, we have expressed 
the matrix elements in a basis of generalized seniority 
s (s counts the number of unpaired nucleons in the nu- 
cleus); |0+) = £,a a |a), ; |0+) = E.ftl*)/- Thc J=° 
terms provide essentially all the contribution to M( 0!/ ) 
that is diagonal in s. The canceling parts, J>0, produce 
almost exclusively cross terms with As=+4. The ma- 
trix elements f{s\M^ 0v ^\s)i are roughly proportional to 
(Smai-s), averaged in parent and grand daughter, while 
the cross terms f{s + A\M^°"' > \s) i are constant -in both 
cases scaled by the larger oscillator quantum number in 
each valence space-. The two body matrix elements of 
the operator M<-°^(J = 0) are almost identical to those 
of the isovector pairing of the nuclear effective interac- 
tion, that is why it acts as a "pair counter". At present 
we cannot offer a similarly simple explanation for the be- 
havior of the J>0 terms. Obviously, when the initial and 
final states have seniority zero, the s=0 contribution is 
maximized and the canceling terms are null, hence, M (0l/) 
becomes maximal. 
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FIG. 2: (color online) The neutrinoless double beta decay 
NME's, denned in equation [2] as a function of the maximum 
seniority of the wave functions. 



These results highlight the role of the seniority struc- 
ture of the nuclear wave functions in the build-up of the 
Qv NME's, and we shall examine this issue for the com- 
peting theoretical approaches. In the first place, we have 
plotted the results of the ISM calculations of the NME's 
as a function of the seniority in Fig. [51 The values with 
maximum seniority provide the exact ISM results in the 
corresponding valence spaces. Two aspects are worth 
to underline; a) the strong reduction of the NME as 
the maximum allowed seniority increases (up to a fac- 
tor five); and b) the fact that, at s<4, the NME's of the 
A=76, 82, 128, and 130 decays miss convergence by fac- 
tors 2-3. On the contrary, in the A=48, A=124, and 136 
cases the convergence at s<4 is much better. The reason 
why these decays behave differently is very illuminating; 

Sn has only neutrons in the valence space, hence, its 
wave function is dominated by low seniority components 
and its NME at s<4 is quite close to the exact result; 
in the A=136 decay, the s<4 calculation for 136 Xe is ex- 
act, therefore, at s<4, the NME is also close to the exact 
one; finally, in the A=48 decay the s>4 components are 
negligible both in doubly magic 48 Ca and in 48 Ti. 




A=76 82 
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FIG. 3: (color online) The neutrinoless double beta decay 
NME's; comparison of ISM and QRPA calculations. Tu07; 
QRPA results from ref. Jy07; QRPA results from 

ref. [13]. ISM s<4 and ISM; present work. 



We can now proceed to compare in detail the "state 



of the art" ISM and QRPA [2l|, |22| NME's in Fig. d 
The QRPA results for 124 Sn are not yet available. The 
range of QRPA values shown in the figure is that given 
by the authors, and derives from the different choices of 
g pp and gA, as well as from their use or not of a renor- 
malized version of the QRPA. The larger values corre- 
spond to g^=1.25 and should therefore be preferred in 
the comparison with our predictions. Both the QRPA 
and the ISM calculations include the higher order correc- 
tions from ref. [loj . For a proper comparison, thc TU07 
NME's have been increased by 10% due to their differ- 
ent choice of rg. In all the calculations the short range 
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correlations are modeled by the same Jastrow factor. 

Several interesting conclusions stem from this figure. 
First, the fact that the different QRPA calculations are 
now compatible. In addition, they produce NME's that 
are strikingly close to the ISM ones calculated at the 
truncation level s<4. In the A=136 decay, in which the 
s<4 truncation is a good approximation to the full result, 
the QRPA values and the ISM ones do a gree (this seems 
to be also the case for the A=124 decay [24|). This sug- 
gests that, somehow, s<4 is the implicit truncation level 
of the QRPA. In the QRPA calculations in a spherical 
basis that we are discussing, the ground states of parent 
and grand daughter, calculated in the BCS approxima- 
tion, have generalized seniority zero. The RPA ground 
state correlations of multipole character (quadrupole, oc- 
tupole, etc.), bring components with s>4 into these wave 
functions. But, for the RPA approximation to remain 
valid, their amplitudes should decrease with s. Indeed, 
in our ISM s<4 results, the percentage of s=0 compo- 
nents is always larger than 70%, a figure compatible with 
a QRPA description. However, in the full calculation 
for the A=76, A=82, A=128, and A=130 decays, this 
percentage can be as low as 25% (actually, in 76 Se, the 
s=4 components almost double the percentage of the s=0 
ones). In these cases, the QRPA is bound to overestimate 
the amount of s=0 components and, consequently, the 
value of the NME's. In a sense, the QRPA can be said to 
be a "low seniority approximation" , roughly equivalent to 
the s<4 ISM truncations, that overestimate the NME's 
when the nuclei that participate in the decay are strongly 
correlated by the multipole part of the effective nuclear 
interaction. The extent of the overestimation depends on 
the degree of validity of the low seniority approximation 
in each decaying pair. 

The values of M 1 - ^ predicted by the present ISM cal- 
culation for the A=76, A=82, A=128, and A=130 de- 
cays, are smaller that the QRPA (central) ones by fac- 
tors 1.5-2. Therefore, for a given value of the effective 
neutrino mass, the predicted ISM half-lives of the Oi//3/3 
decays are 2-4 times longer than the QRPA ones. Equiv- 
alently, for a given lower bound on the half-life, the ISM 
NME's produce upper bounds on the effective neutrino 
mass that are larger than those of the QRPA by factors 
1.5-2. For instance, a bound on Ti ( 76 Ge -»■ 76 Se) of 10 25 
y. results in an effective neutrino mass of 1.05 eV with 
the ISM NME, and 500 meV with the QRPA one. The 
same bound for the half-life of the 130 Te — ► 130 Xe decay 
would lead to bounds on the neutrino mass of 410 meV 
and 270 meV respectively. 

In summary, we have analyzed the 0v(3(3 NME's in 
terms of the J 77 of the decaying neutron pair. We have 
found that in the seniority zero limit the decays are 
strongly favored. When the non zero seniority compo- 
nents of the wave functions, originated by the multi- 
pole terms of the nuclear effective interaction, are prop- 
erly taken into account, the matrix elements are dras- 



tically reduced. In particular, when the multipole cor- 
relations are large, the low seniority truncations, s<4, 
similar to those implicitly present in the spherical QRPA 
approaches based in a BCS treatment of the pairing in- 
teraction, are shown to overestimate the NME's. Hence, 
we surmise that, when the QRPA and ISM results do not 
agree, the true NME's should be much closer to the ISM 
predictions than to the QRPA ones. 
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